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We present the solution of the Poincare-covariant Faddeev equation for the A(1232) and f2(1672) 
baryons. The covariant structure of the corresponding baryon amplitudes and their decomposition in 
terms of internal spin and orbital angular momentum is explicitly derived. The interaction kernel is 
truncated to a rainbow-ladder dressed-gluon exchange such that chiral symmetry and its dynamical 
breaking are correctly implemented. The resulting physical masses agree reasonably with experiment 
and their evolution with the pion mass compares favorably with lattice calculations. Evidence for 
the non-sphericity of the A— resonance is discussed as well. 
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I. INTRODUCTION 

Quantum Chromodynamics (QCD) was built upon the 
understanding of hadrons via the quark model and the 
discovery of color. However, despite the success of QCD 
predictions at high energies, which led to the consensus 
that QCD provides a correct description of strong inter- 
actions, the calculation of hadron properties developed 
at a slower pace. This is due to the fundamental non- 
perturbative character of bound-state phenomena, and 
it is further complicated by the low-energy enhancement 
of strong interactions. 

Hadrons, as bound states of quarks and gluons, are 
described in continuum QCD by (generalized) Bethe- 
Salpeter equations (BSEs) [1-5] which rely upon the 
Green functions of the theory. This, in principle, neces- 
sitates a solution of the infinite tower of coupled Dyson- 
Schwinger equations (DSEs) [6-8]. Any feasible numeri- 
cal procedure thus requires a symmetry-preserving trun- 
cation, but once such a truncation is performed, the ap- 
proach provides a fully consistent quantum-field theoret- 
ical framework for the study of hadron properties. 

The aforementioned approach has a longstanding and 
successful history in the investigation of meson proper- 
ties; see, e.g., [9-11] and references therein. For baryons, 
however, the relativistic bound-state description is con- 
siderably more involved and computationally demand- 
ing. The reason is that the presence of a third quark 
enlarges the momentum phase space, and the relativis- 
tic spin structure [12, 13] of the three-quark bound-state 
amplitude is much more complicated as well. A success- 
ful simplification of the problem starts from a covari- 
ant Faddeev equation [3, 4] but immediately reduces its 
complexity by treating the two-body scattering matrix in 
a separable expansion: by considering quark-quark cor- 
relations, called diquarks, as dominant, the three-quark 
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equation is simplified to a quark-diquark BSE [14-16]. 
For a collection of recent results using the quark-diquark 
approach see, e.g., Refs. [17-23]. 

A step forward was taken in Refs. [13, 24-26] where, for 
the first time, the full three-body Poincare-covariant Fad- 
deev equation for the nucleon was solved. Here, the in- 
teraction kernel is truncated to a rainbow-ladder dressed- 
gluon exchange between two of the quarks, and irre- 
ducible three-quark contributions are neglected. The nu- 
cleon mass reported in these works is comparable to lat- 
tice data and, surprisingly, also very close to the quark- 
diquark result. The calculation has been recently im- 
proved in Ref. [27]. 

In the family of baryon resonances, the A(1232) iso- 
quadruplet plays a special role, owing to its high produc- 
tion cross section which makes it one of the best studied 
resonances. Moreover, it is the lightest baryon resonance 
and possesses the same valence-quark content as the nu- 
cleon, and thus constitutes an excellent system to under- 
stand excitations in QCD. It is also the lightest known 
particle with spin 3/2- At present, due to numerical lim- 
itations, such covariant three-body bound-state calcula- 
tions are performed only for baryons with equal-mass va- 
lence quarks. This is the case for the A quadruplet as- 
suming isospin symmetry. Another spin-3/2 ground-state 
baryon is the 17(1672): it is a pure strange state and 
therefore allows to study the quark-mass dependence of 
A properties. 

In the present work we extend the techniques described 
in [13, 25, 27] to the case of the A and fi resonances in 
view of computing their masses and bound-state ampli- 
tudes. We analyze the importance of the different partial- 
wave contributions and thereby demonstrate that the A 
and VL baryons are not pure s-wave states. This repre- 
sents a first step towards the understanding of baryon 
deformation from sphericity. A more complete descrip- 
tion of baryon structure will require to study the con- 
tributions of different internal spin and angular momen- 
tum structures to physical observables such as, e.g., elec- 
tromagnetic form factors or the electromagnetic 'yN IS.- 
transition. 

The paper is organized as follows. In Section II 



2 



we review the Poincare-covariant Faddeev approach to 
baryons and the rainbow-ladder truncation of the inter- 
action kernel, and we describe the main features of the 
covariant decomposition of the Faddeev amplitude for 
spin-3/2 baryons. In Section III we present our results for 
the A and Q masses and compare with lattice data and 
experimental values. We also discuss the internal angu- 
lar momentum composition of the amplitudes. Techni- 
cal details about the partial-wave decomposition and the 
numerical implementation are described in Appendices 
A and B, respectively. Our calculations are performed 
in Euclidean momentum space and Landau gauge QCD; 
details about the conventions used in this paper can be 
found, e.g., in Ref. [27]. 



II. FADDEEV EQUATION AND 
RAINBOW-LADDER TRUNCATION 



Baryons are described in QCD by the three-quark (six- 
point) Green function which characterizes the state both 
on- and off-shell. On the baryon mass shell the Green 
function has a pole, and a Laurent expansion around this 
pole allows to derive a relativistic three-body bound-state 
equation (see Ref. [5] for a pedagogical discussion). That 
equation reads 



the physical normalization condition 



a=l 



(a) 
(2) 



(1) 



where 5" represents the baryon's bound-state amplitude 

(3) 

and Gq is the product of three dressed quark propaga- 
tors. The three-body kernel includes all possible cor- 
relations among the three quarks: it comprises a three- 
body irreducible contribution K^^^ as well as the sum of 
the three two-body irreducible interactions with a specta- 
tor quark a, -^^[2)- The bound-state amplitude ^' satisfies 



d 



G. 



(3) 



(2) 



where ^ is the charge-conjugated amplitude and P the 
total baryon momentum. 

The structure of the bound-state amplitude depends 
on the baryon of interest. It is the product of color, fla- 
vor and spin parts. For spin-i/2 baryons, the spin part 
is a rank-four Dirac tensor, with three Dirac indices rep- 
resenting the valence quarks and the fourth index the 
bound state. A spin-3/2 particle is described by a Rarita- 
Schwinger field so that in this case the spin part of the 
bound-state amplitude is a mixed tensor with four Dirac 
indices and one Lorentz index. Moreover, it depends on 
the three quark momenta pi , p2 and pa which can be con- 
veniently reexpressed in terms of the total momentum P 
and two relative momenta p and q: 



P = (1 - C)P3 - (Pd , Pi 
P2 -Pi 

q = — ^ — , P2 = q 

P =Pl+P2+P3, P3 



P 



p 

2^^ 

p + CP, 



p. 



p, 



(3) 



where Pd — Pi + P2- The total momentum is constrained 
by P^ = —AP, with M being the baryon mass. C is 
a free momentum partitioning parameter which we will 
choose to be C = V^- This choice maximizes the accessi- 
ble bound-state mass range with respect to the analytic 
structure of the quark-propagator [24] and also allows for 
a tremendous simplification of the bound-state equation's 
solution method (see Appendix B). 

The success of the quark-diquark approach to baryon 
properties supports the idea that quark-quark correla- 
tions dominate the binding of baryons. Consequently, 
we neglect the three-body irreducible contribution Ki[^ . 
Then, considering the above discussion about the bound- 
state amplitude's tensorial structure and kinematics, one 
arrives at the covariant Faddeev equation for a spin-3/2 
baryon (see Fig. 1): 



Kp^sip^i^P) 



Kp0.,,. {k)Sp.p., {k2)Syy, {h)Kp"Y's{P^'^ P) 
+K^Yo.o.' ik)Syy, (fc3)^a'a" {~kl)K" sip''^^ . 9^'^ : P) 
+K^c.'PP'{k)S„,^,,{k^)Sp,p„ik2)K"P"jsiP^'\<l^'^^P) 



(4) 



where we have already restricted ourselves to the case 
where the two-quark kernel depends only on the exchange 
(i.e., gluon) momentum k. The quark propagators S de- 
pend on the internal quark momenta ki = pi — k and 



ki — Pi + k. The internal relative momenta are 

= p k, p^^-* ^ p — k, p^^-* = p, 
g(i)=g-fc/2, q^^^^q-k/2, q^^^^q + k. 



The solution of the Faddeev equation is greatly simpli- 
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FIG. 1. Diagrammatical representation of the Faddeev equation in rainbow-ladder truncation. 



fied if one expresses the Faddeev amplitude in terms of a 
basis T'^f3^s''{p,q,P), 



where the expansion coefficients fi depend on the five 
Lorentz-invariant variables 



zi 



p-P, Z2^q-P. (7) 



Here, a hat denotes a normalized four-vector and the 
subscript T a transverse projection with respect to the 
total momentum P. 

In principle one is free to choose any possible basis r*. 
For the sake of physical interpretation it is convenient 
to perform a partial-wave decomposition and classify the 
basis elements with respect to their quark-spin and rel- 
ative orbital angular momentum content in the baryon's 
rest frame. The explicit derivation of such a basis is quite 
involved [28] and presented in Appendix A. The most sig- 
nificant aspects of this construction are: 

• It is independent of any approximation in the Fad- 
deev equation. 

• Only Poincare covariance as well as parity invari- 
ance are needed to construct the basis. 

• For positive-parity, positive-energy (particle) spin- 
3/2 baryons the basis consists of 128 elements. 

• The basis includes all possible spin and orbital an- 
gular momentum values, namely, s — 1/2 , 3/2 and 
£ = 0,1,2,3. 

Once the Poincare-covariant structure of the spin-3/2 
baryon amplitude is determined, the elements in the Fad- 
deev equation that remain to be specified are the two- 
quark interaction kernel and the quark propagator. They 
are related by the axial-vector Ward-Takahashi identity 
(AXWTI) which, in pseudo-scalar meson studies, guar- 
antees the correct implementation of chiral symmetry 
and its spontaneous breaking. In particular, it implies 
that pious are massless in the chiral limit [29-35]. The 
simplest qq interaction kernel that satisfies the AXWTI 
is a dressed-gluon ladder exchange: 



ICaa'PP'ik) 



where Z2 is the quark renormalization constant, T^" — 
(5'"' _ k^'-k'^ is a transverse projector with respect to the 
gluon momentum fc, and a{k'^) is an effective interaction 
that defines the model input. 

The quark propagator satisfies the quark Dyson- 
Schwinger equation whose interaction kernel includes 
the dressed gluon propagator as well as one bare and 
one dressed quark-gluon vertex. Consistency with the 
AXWTI requires to use the same interaction (8) in the 
quark DSE: 

^ap(P) = ^2 (i^ + ™o)a/3 + / ^aa'f5'fs{k) Sa'P'{q) , (9) 



where toq is the bare current-quark mass that enters the 
equation as an input, and k = q — p. Tracing the color 
structure yields a prefactor 4/3 in front of the DSE inte- 
gral in (9) and 2/3 for the integral in the Faddeev equa- 
tion (4). Eqs. (8-9) define the rainbow-ladder trunca- 
tion where the dressed quark-gluon vertex is truncated 
to its vector part 7^ , and the combined non-perturbative 
dressing of the gluon propagator and quark-gluon vertex 
is absorbed in the effective interaction a{k^) which has 
to be modeled. 

For the effective coupling we adopt the ansatz intro- 
duced in Ref. [36] which has been successfully applied 
in many hadron studies. It reproduces the logarithmic 
behaviour of QCD's one-loop running coupling in the ul- 
traviolet, and it is strong enough in the infrared to enable 
dynamical chiral symmetry breaking and thereby gener- 
ates a dynamical constituent-quark mass: 

a(fc2) ^ nrj' (^'f^jV^^'^' /^') + auv (k^) ■ (10) 

The infrared part is characterized by one dimensionless 
parameter rj and an energy scale A. The 'ultraviolet term' 
reads 



auv(fc^) = 



27r7,.(l-e-^-V^?) 
ln[e2-l + (l + fcVA|cD)']' 



(11) 



k la 



■T/3/3' 



(8) 



where 7„j — 12/{llNc — 2Af/) is the anomalous dimen- 
sion of the quark propagator. We use 7,„ ~ 12/25 and 
Aqcd = 0.234 GeV. The scale A* = 1 GeV is only in- 
troduced for technical reasons, its precise value has no 
impact on the results. The infrared energy scale is fixed 
to the value A = 0.72 GeV in order to reproduce the pion 
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Nucleoli 


A 


n 


Physical 


0.94 


1.23 


1.67 


Q-DQ [23] 


0.94 


1.28 


1.77 


Faddcev 


0.94 [27] 


1.26 


1.72 



TABLE I. Nucleon, A and Q masses (in GeV) obtained from 
the Faddeev equation and, for comparison, in the quark- 
diquark (Q-DQ) approach [20, 23]. In both cases, the pa- 
rameter rj in Eq. (10) is set to 77 = 1.8. 



decay constant [36]. Many meson observables, and here 
especially the ones related to ground-state pseudoscalar 
and vector mesons, turn out to be almost insensitive to 
the infrared width parameter 77 around a central value 
r] = 1.8 [11], and a similar observation holds for nucleon 
properties [23, 27]. 




0.2 0.3 
mn^ (GeV^) 



III. RESULTS 

Having defined the input of the Faddeev equation, 
Eqs. (4) and (9) can now be solved consistently. The de- 
tails of the calculation are given in Appendix B. From the 
solution of the Faddeev equation we obtain the bound- 
state mass and amplitudes, and the Faddeev amplitudes 
thereby provide information about the internal structure 
of the baryon. 

For the hadron mass calculation, the current-quark 
masses must be fixed as well. Upon setting the value 
for the infrared scale A, the u/d mass is adjusted to re- 
produce the physical pion mass. Since there is no purely 
strange pseudoscalar meson, the s-quark mass is fixed to 
reproduce the kaon mass which would correspond to a 
fictitious pseudoscalar ss state iriss — 690 MeV. 

To analyze the sensitivity of our results to the param- 
eters of the effective interaction, we have repeated our 
calculations for different values of 77 in the range 1.6 — 2.0. 
For the central value rj = 1.8, the resulting A mass at the 
physical pion mass, and the fl mass at ttItt = 690 MeV, 
are shown in Table I. They are in good agreement with 
the corresponding experimental values. For comparison, 
we have also included the respective mass results in the 
quark-diquark approach [20, 23]. While the same value 
for the nucleon mass, Mjv = 0.94 GeV, is obtained both 
in the quark-diquark approach and through the Fad- 
deev equation, the quark-diquark results for the spin-3/2 
baryon masses are larger. This solidifies the observation 
that a truncation of the Faddeev kernel to (ground-state) 
scalar and axial-vector diquark channels is sufficient to 
reproduce nucleon properties, whereas more structure is 
needed to describe spin-3/2 baryons. Nevertheless, the 
results clearly demonstrate that diquarks capture an im- 
portant part of the internal structure of baryons. 

It is interesting to study the relevance of interaction 
terms beyond the rainbow-ladder truncation. Meson 
studies show that pionic corrections (pseudoscalar-meson 
exchange among quarks, so-called resonant corrections) 



FIG. 2. Evolution of Ma with compared to lattice data 
[37-40]. Stars denote the experimental values of A and n. Er- 
ror bars in the Faddeev calculation correspond to a variation 
in the interaction width in the range rj £ 1.6 — 2.0 



are attractive, while non-resonant contributions have a 
repulsive effect [41-44]. Both effects seem to cancel at 
the non-perturbative level, so that the rainbow-ladder 
kernel gives an accurate description of the qq interaction 
in mesons. Moreover, pionic corrections are expected to 
be suppressed at high quark masses. 

The evolution of the A-mass with the squared pion 
mass is shown in Fig. 2 and compared to lattice results. 
Here, the pion mass is obtained by solving the corre- 
sponding Bethe-Salpeter equation using the same inter- 
action kernel upon varying the current-quark mass. We 
find a good agreement of the Faddeev-calculated masses 
with lattice data, especially at higher pion masses. If 
the analysis of meson results for effects beyond rainbow- 
ladder can be extended to the gg-interaction in baryons, 
this would indicate that a near cancelation of reso- 
nant and non-resonant contributions prevails up to the 
strange-quark mass and beyond. 

Ultimately, such issues can only be clarified by over- 
coming the rainbow-ladder truncation as well as the 
restriction to quark-quark correlations in the Faddeev 
equation. We note, however, that the current-mass evo- 
lution of the A mass from the Faddeev equation is in 
qualitative agreement with that of the vector-meson mass 
obtained from its Bethe-Salpeter equation, cf. Fig. 2 in 
Ref. [23]. Using the current-mass independent interac- 
tion defined in Eq. (10), both resulting curves overes- 
timate the experimental values of Mq and m$ at the 
strange-quark mass approximately by the same amount 
(^ 3%). In view of this, it could be misleading to at- 
tribute such discrepancies solely to corrections beyond 
rainbow-ladder without understanding the sensitivity of 
our calculations to the effective interaction a{k'^). It 
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would be desirable to compare results using different, 
lattice-inspired, model interactions; see also Ref. [45]. 
Such a study will be published elsewhere. 

There is another missing feature in the current-mass 
evolution of the A mass. The A resonance can decay 
into a nucleon via the emission of a pion. Such a de- 
cay would manifest itself in a non-analytical behaviour 
of the A mass as a function of the pion mass when the 
decay channel opens, i.e. for ttItt ~ Ma—Mn ~ 300 MeV 
[46, 47]. Our truncation scheme does not provide a mech- 
anism for the A to decay and thus we should obtain 
a higher A mass, shifted approximately 100 MeV up- 
wards, corresponding to the A decay width. Again, an 
analogous observation holds for the p— meson obtained 
through the Bethe-Salpeter equation [23, 36]. Therefore, 
we conclude that the model provides too much binding 
in the light pion-mass region where the decay channel 
would be opened in a full calculation. 

A first hint on the internal structure of the A-baryon 
comes from the analysis of the relative importance of the 
different quark-spin and relative angular-momentum con- 
tributions. In Fig. 3 we plot the dominant amplitudes in 
each partial- wave sector (see Appendix A). It is clear that 
the s = 3/2, £ = [i.e., s-wave) sector is dominant. How- 
ever, p- and d-wave sectors are much larger in number of 
basis elements, thus even small values of individual basis 
elements can add up to a non-negligible contribution of a 
given sector and thereby cause deviations from sphericity. 
We obtain very similar results for the Jl-baryon. In any 
case, a thorough understanding of the A and internal 
structure requires the study of the different partial- wave 
sectors in physical observables, e.g., form factors. 
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FIG. 3. Zeroth Chebyshev moments of the coefficients fi in 
the basis expansion of Tables II and III. We plot the coeffi- 
cients dominant in each (i, ^)-sectors as a function of p and q. 
The normalization is chosen such that fi{p^ = 0,q^ = 0) = 1. 
The respective basis elements t^^' are from left to right, top 
to bottom: t+/, t+\ t^^^, t;^^\ t+^, 



IV. CONCLUSIONS 



We presented a fully Poincare-covariant solution of the 
Faddeev equation for the A- and f2-baryons. The re- 
quired quark propagators are obtained by consistently 
solving the corresponding Dyson-Schwinger equations. 
For this we employed a rainbow-ladder truncation of 
the interaction kernel, dressed by an effective interaction 
which effectively depends on a single parameter and is 
fixed to reproduce the pion decay constant. 

We obtain masses that are in good agreement with ex- 
perimental results and their evolution with the pion mass 
compares favorably with lattice results. We concluded, 
however, that the effective interaction we used provides 
too much binding in the low pion-mass region. 

The contribution of the different quark-spin and rel- 
ative orbital angular momentum sectors was studied by 
naively comparing the magnitude of the corresponding 
dominant amplitudes. This shows an s-wave dominance, 
however with a non-negligible p- and d-wave contribu- 
tion. A better understanding of the internal structure 
would come from the study of A and fl form factors 
in the present framework. Work in this direction is in 
progress. 
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Appendix A: Covariant decomposition of the 
A— baryon amplitude 

The spin part of the Faddeev amplitude for spin-3/2 
particles is a mixed tensor with four Dirac indices and one 
Lorentz index. The positive-parity and positive-energy 
subspace can be expanded in a basis with 128 elements 
T^o . ^(p, g, P), cf. Eq. (6), which can be classified with 
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3S?4 + 2Vfi 


s 
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J_ <T,1 


3/2 
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S?i + S5i+2S|i-|(Vl3 + 2V?4) 
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2 


S?2-2SSi-|(Vf3-2Vl4) 
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•^^13 ' \'-'43 ^64:1 ' 3 V 11 T V 

S?4-2(S^3-S54) + i(V?i-2Vl2) 
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3/2 
3/2 


2 
3 


V5(S?i+3S5i-Vl3) 

4S?2+5(S52 + S|2)+V'i4-V?3) 


3/2 


3 


Sx3 ~t~ 5 S33 Vii 


3/2 


3 


Si4 + 5 S44 — Vii 


s 


I 


J_ <T,1 


3/2 
3/2 
3/2 
3/2 


3 
3 
3 
3 


S?i + 2 (Si2 + Sli + S5i) - 1 (V'i3 + 2 Vf4) 

%2 - SI2 - 1 (V?3 + Vl4) 
S?3+S53 + 2S|3+|(V'n+2V?2) 
S?4 + S|4 + 2S54+|(V?l-2Vl2) 



TABLE II. Orthonormal Dirac basis rf-^ for 5 — 3/2 and 
k — 1. We omit Dirac and Lorentz indices as well as the label 
a = ± for better readability. The dominant covariants in each 
(5, i')-sector which are plotted in Fig. 3 are highlighted; the 
relevant values of a and k are given in the caption of Fig. 3. 



respect to their quark-spin and relative orbital angular 
momentum in the baryon's rest frame. The resulting 
basis is collected in Tables II and III, and in the following 
we will sketch its derivation. 

To begin with, it is convenient to express the momenta 
{p, (7, P} that enter the basis elements through orthogo- 
nal unit vectors {r, s, P} which satisfy = = = 1 
and r ■ s — r ■ P — s ■ P — Q. This is achieved by taking 
the component of p transverse to P, and by projecting q 
onto the direction transverse to P and r: 



Pt 



P - ZiP 



qr 



Z2P 



with the Lorentz-invariant momentum variables zi, Z2 
and zq from Eq. (7). In the baryon's rest frame and a 
suitable choice of coordinates, r, s and P become the 
Euclidean unit vectors 63, 62 and 64. 
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J_ a,l 


(T,l 
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VJ3 + VU 
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1/2 
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* 11 
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Sii 


Vii 
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CT.l 


V3</ 


1/2 
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or QS 
313 ^ 'J14 


VLs - v?4 


1/2 
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S13 + Si4 


v?3 + v;4 


1/2 
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Sii + 2 S12 


V'n + 2V?2 


1/2 


2 


OS n or 
Bii — Z ai2 


V?i-2Vl2 



TABLE III. Orthonormal Dirac basis t'^''' for 5 = V2 and 
fc = 1. We omit Dirac/Lorentz indices and the label a. The 
dominant covariants in each (5,^)-sector are highlighted. 



Next, we define the basic Dirac structures 

r.=i...4 - {1, f, ^} (Ai) 

which we use to construct a (still linearly dependent) 
basis for a rank-four Dirac tensor: 



( 



1 (g) 1 

7^ ® 7^ 



\ 



{T,(g)Tj)n''{P). (A2) 



V 7^7' ® 7^7' / 



Here, 7^ is the 7— matrix transverse to P, and we used 
f7±(P) = A±(P)75C (g) A+(P), where C = -f'^-f is the 
charge conjugation-matrix and A*(P) — (1 ± /2 are 
the positive- and negative-energy projectors. The tensor 
products are understood as 



{h® f2){9l®92) = (/igi)«)(/2 52) 



(A3) 



Denoting the elements in Eq. ( A2) generically by M^^ , 
with M S {S,P,V, A} and «,j = 1...4, we can con- 
struct the building blocks of a Dirac-Lorentz basis for 
the A— baryon in the following way: 







( 7^ 7^ (g) 1 ^ 




[Mir 




r''7^ (g 1 








s'^75 g) 1 y 





where f'^'^ is the Rarita-Schwinger projector for positive- 
energy particles: 



'\P)=K+{P) (t^'' 



1 



(A5) 
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The set (A4) contains 384 elements, but one can show 
that only 128 of them are linearly independent. We found 
the following choice of linearly independent elements con- 
venient: 







'^43 , ^41 


pr -ps 




'^32 1 ^34 ^ 






nr pr 
-t^43 , '^41 






ps ps 
L J^^32 1 J^^34 J 



(A6) 
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, 
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with j = 1 ... 4, and we have omitted the index a for 
better readability. The set (A6) includes the element 

(A7) 



which is familiar from the quark-diquark model and rep- 
resents the dominant s— wave structure in the A— baryon 
amplitude. 

Let us now turn to the partial-wave analysis of these 
basis elements. In the baryon rest frame, the quark total 
spin and relative angular momentum operators read 



5-2 ^ 1(1 (g) 1 (g) 1) 



7^2 



(</) 



2 L/„\ ■ L 



crj. (g) 1 -I- perm.) 



where al^'^ 



r2 



^(p) 

'i [7t,7t] and 



(9) 



(A8) 



2pVp+pnpVp)V^p-p2Ap, 



(<]) 



= 2q.Vq + g'(q- Vq)Vjj-q2Aq, (A9) 

hp) ■ hi) = • Vp)Vq - (p • q)(Vp • ^q) , 

and p, q arc the spatial parts of pt and qt, respectively. 

It is useful to realize that the basis elements containing 
{S, V} and {P, A}, which differ by a factor 7^ g) 7^, and 
those with a different value for the index cr = ±, do not 
mix under the action of or and thus can be analyzed 
independently. Moreover, from Eqs. (AS) and (A9) one 
infers that the set (A6) can be further subdivided into 
four subsets which, due to their different momentum de- 
pendence, again do not mix under or L^: 

1 „2 2 2 2 . c9 qr OS or os os y-r tts 
i^, I , i> , I b . Dii, 0]^3, O44, O34, O44, O32, \ 12, V 44 

r„ ^3 3 . OS or OS OS or or y-r ys 

I t>, I b, I t> . O42, i->i4, 'J13, O42, O44, O32, \ 14, V 43 

■ '-'13' ^11' '-'12, ''33' '^34' '-'43' ''ll' ^12 

„ ^2 3 . C9 or OS os or or yr ys 

i>, I b, t> . O44, D12, Oil, ^44' '-'34, '^43' *12i * 11 



The operator S*^ is independent of the momentum con- 
tent of the basis elements. Therefore it is sufficient to 
find the eigenstates at fixed values for the unit vectors 
r and s. At this point the problem can be easily imple- 
mented and solved using a symbolic programming lan- 
guage. Similarly, the orbital angular-momentum decom- 
position can be done straightforwardly on a computer, 
but it is nevertheless instructive to study some simple 
examples. The i = elements can be found immediately: 
they are the momentum-independent elements in (A6), 
i.e., [Sfi]'^ and [Pfi]'^. The remaining basis elements can 
be written as contractions of 



r", s", r"s'^, r^rf^, s°'s^^, r'^r^s^ 



(AlO) 



with appropriate momentum-independent Dirac-Lorentz 
tensors. For the £ = 1 elements it is sufficient to consider 
the first three elements in the list above. Applying the 
orbital angular-momentum operator yields 



L2s" = 2 s" , 

and thus the £ = I cigcnfunctions are given by: 
r , s , r s — s . 



(All) 



(A12) 



Their contraction with the corresponding Dirac-Lorentz 
structures yields the £ = 1 eigentensors at the level of 
the basis elements. For other £ values the calculation 
proceeds along the same lines but is more involved. 

The above analysis focused on the subset {S, V}, which 
we denote by the index fc = 1. The final result of the 
partial-wave decomposition for fc = 1 is presented in Ta- 
bles II and III. Analogous results hold for the set {P, A}, 
denoted by k=2: the respective basis elements are ob- 
tained by exchanging S — J' P, V — >■ A and adding an 
extra minus sign to the elements Py . For the compi- 
lation in Tables (II-III), the original label i — 1 ... 128 
that appears in Eq. (6) was dissolved into the indices 
i = 1 . . . 8, j = 1 . . . 4, fc = 1, 2, and CT = ±. 

The resulting basis satisfies the following orthonormal- 
ity relation: 



:Tr. 



(A13) 



\ ij i'j' j si' ij ifiaSjli'j' i afifS 
= Sii'Sjj'Skk'Sacr' , 

where the conjugation of the basis elements is defined as 



11' ^^'a'S'yi^P^ ^J'^ ^I'6 ■ 



(A14) 



Here, the left column symbolically indicates the different 
momentum dependencies of the basis elements in powers 
of r and s. The partial-wave analysis is now considerably 
simplified since we only need to look for eigenfunctions 
of S"^ and within the above subsets. 



Appendix B: Numerical details 

The numerical techniques used in this work are an ex- 
tension of those explained in Ref. [27] towards the case 
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FIG. 4. First four Chebyshev moments in the variables zq, 
Zi and 22 for the dominant amplitude [iSfj^]^ at — 0. Note 
that the zeroth order is the same for all three figures; in the 
lower two panels it is not shown for clarity. 



of the A-baryon. We summarize the main ideas in this 
appendix. 

The Faddeev amplitude is the product of color, flavor 
and spin parts. As required by the Pauli principle, that 
product must be antisymmetric under the exchange of 
any two of the three quarks. Since a baryon is a color 
singlet, the color part of the Faddeev amplitude is always 
antisymmetric, and thus the product of flavor and spin 
parts must be symmetric. This symmetry requirement 
allows to relate the three diagrams in Fig. 1. 

Compared to the nucleon, the situation is simpler for 
the A- and Sl-baryons. Since the A carries isospin 3/2, its 
flavor part is completely symmetric and hence also the 
spin part must be symmetric. The SI, on the other hand, 
is a pure s-quark state and thus automatically symmetric 
in flavor space. Using this, and starting from Eq. (4), it 



is not difflcult to show that the first two terms in the 
equation can be reexpressed in terms of the third one, 
evaluated at different momenta: 



-.(3) 



(Bl) 



''"^^ 'To" n" P) 



where the permuted relative momenta are given by 



p = -q- 



9' =-! 



3p 
4 
3p 

4 



(B2) 



In terms of the amplitude decomposition (6) , the Fad- 
deev equation is written as 



/,(p,g,P)-^/('^'(p,g,P) 



(B3) 



where we used the notation 

(p, q, P) = / (p, q, P; k) (pf"^) , qf-^) , P) 



4,'^t\p,q,P)^Q[-\p,q,P)h{p,q,P): 



(B4) 



with p^""^ and g^"-* being the internal relative momenta 
from Eq. (5). The functions (j)f'^ are the coefficients of 
the so-called wave functions ^^^^ = S{pb)S{pc)'^ in the 
basis expansion (6), where {a, 6, c} is a symmetric per- 
mutation of {1,2,3}, and ICij and Qij are the matrix 
representations of the kernel and propagator operators 
in that basis. Eq. (Bl) can be written in terms of the 
coefficients fi{p,q,P) as 



(B5) 



Mp,q,P) = f^'\p,q,P) 

+ Hl^ff\p',q',P) 
+ ^ff\p",q",P), 

where the matrices H' and H" are defined by 

Hl^ ip,q,P) = l f^„,, b, g, P) T^^^s " iP' ^I'^P). 

H^ij,, q,P)^\ f^„,, '^b, q, P) r^^ps'iP"^ l" ^ P) ■ 

(B6) 

In rewriting the equation we have used the orthogonal- 
ity (A13) of the basis. 

The algorithm for solving the Faddeev equation begins 
by introducing a fictitious eigenvalue \{M) in front of the 
integral in (B4). The integral equation is then solved by 
iteration: start with a guess for the bound-state mass M 
and the initial amplitudes fi] compute the integral (B4) 
for the case a = 3; apply the permutation (B5); and 
repeat until the eigenvalue X{M) has converged. If the 
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converged value is A(M) = 1, then M is the correct mass 
for the bound state; otherwise one has to change the 
guess for M and repeat the procedure. 

For the presented solutions of the Faddeev equation we 
expanded the angular dependence of the Faddeev ampli- 
tudes in terms of Chebyshev polynomials. This is con- 
venient because only a small number of Chebyshev mo- 
ments are needed, as shown in Fig. 4. 

Solving the Faddeev equation in the form (B4-B5) is 
a great simplification with respect to the original formu- 
lation (4) since the kernel matrix /C*^^^ is comparatively 
simpler than the other two permuted kernels. Note that 
in this case the external and internal relative momenta p 
and p^'^^ are the same, cf. Eq. (5). If one uses a moder- 



ate number of integration points, and taking into account 
the fact that the kernel matrices are typically sparse, then 
the kernel matrix can be computed and stored in advance 
which notably reduces the computing time. However, if 
the number of integration points is increased, memory 
limitations become an issue. For the presented results 
we have used the following number of integration points: 
20 for {p^,9^,fc^}, 8 for zq and 4 for {zi,Z2,z,y}. For 
this configuration, and without storing the kernel matrix 
in advance (although it would be feasible), each itera- 
tion requires approximately four CPU-hours on a 16-node 
2.66 GHz cluster, and convergence for \{M) is usually 
reached within 15-20 iterations. 
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